Superconductivity from a long-range interaction: 
a crossover between the electron gas and the lattice model 
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We explore how the superconductivity arising from the on-site electron-electron repulsion will 
change when the repulsion is changed to a long-ranged, 1/r-like one by introducing an extended Hub- 
bard model with the repulsion extending to distant (12th) neighbors. With a simplified fluctuation- 
exchange approximation, we have found for the square lattice that (i) as the band filUng becomes 
dilute enough, the charge susceptibility becomes comparable with the spin susceptibility, where p 
and then s pairings become dominant, in agreement with the result for the electron gas by Takada, 
while (ii) the d-wave, which reflects the lattice structure, dominates well away from the half filling. 
All these can be understood in terms of the spin and charge structures along with the shape and 
size of the Fermi surface. 
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While the discovery of high-Tc cuprates has kicked off 
intensive studies of electron mechanisms of superconduc- 
tivity, these have been primarily focused on the pairing 
from short-ranged electron-electron repulsion. This is a 
reasonable assumption for transition metal compounds, 
since the main interaction between the d electrons should 
be short-ranged, as captured in the Hubbard model with 
an on-site repulsion. However, if we go back to the his- 
tory of electron mechanisms of superconductivity, there 
is an important predecessor — the electron gas with the 
long-range Coulomb interaction, one of the most funda- 
mental problems in the condensed-matter physics. The 
problem has a long history, where Kohn and Luttinge ri, 
as early as in the 1960's, have suggested that normal 
states in the electron gas should become unstable in favor 
of a superconducting state, which has later been proved 
when the gas is dilute enough^. There are two important 
factors that discriminate the electron gas from electron 
systems with short-range interactions. One is that the 
gas becomes more strongly interacting (i.e., larger ratio 
of the interaction to the kinetic energy) for the more di- 
lute concentrations. More importantly, the long-range 
interaction can make the charge fluctuations as large as 
the spin fluctuations as contrasted with systems with 
short-range interactions where the spin fluctuations dom- 
inate, so that the fluctuation-mediated interaction should 
change with the range of the interaction. 

A fascinating question then is: what will become of 
the superconductivity in the electron system with short- 
range repulsions when the interaction range is increased 
to approach 1/r. Besides the problem of the interac- 
tion range, there is another essential question: if we 
consider a lattice systems as in the Hubbard model, 
the band filling is a crucial parameter with the half- 
filling being a special point, which controls Mott's metal- 
insulator transition as well as the Fermi surface nesting, 
which in turn dominates the spin fluctuation. Indeed, 
we have a d-wave pairing mediated by spin fluctuations 



in the on-site Hubbard model around the half filling, as 
theoretically suggested with, among other methods, the 
quantum Monte Carlo method^, the fluctuation exchange 
(FLEX) approximation^' and the dynamical cluster 
approximation^. By contrast the electron gas in a contin- 
uous space, only characterized by the electron concentra- 
tion (or, more precisely, r^, the mean electron separation 
measured by the Bohr radius) with a circular Fermi sur- 
face, has no such special fillings. Takada^ has intensively 
studied superconductivity in the three-dimensional elec- 
tron gas, and has concluded that a p-wave pairing should 
occur for Vg > 3.3, which gives way to an s-wave pairing 
when the gas becomes more dilute (r^ > 8.6), where the 
pairing is interpreted to arise from charge fluctuations 
including plasmons^^. 

One step toward the increased range of interaction is 
to consider the extended Hubbard model which takes ac- 
count of the nearest-neighbor repulsion V. The extended 
Hubbard model has been studied with various theoretical 
methodn^^i^^i^'^i^^i^-^^^i^^i^^i^^i^"i^^i^^i^-^i^l For instance, 
Refs l20ll2ll address, with FLEX, the question of how the 
enhanced charge fluctuations in the presence of V should 
affect the pairing symmetry, where the result indicates 
that a triplet pairing as well as singlet ones appear in the 
phase diagram, whose mechanism can be traced back to 
the structure in the charge and spin susceptibilities. 

Given this background, the purpose of the present work 
is to explore what happens to the pairing symmetry as 
we make the interaction more long-ranged and closer to 
the Coulombic interaction. For this purpose we take a 
model where the interaction is extended to distant (12th 
here) neighbors for the square lattice, which is studied 
with a simplified FLEX. We then question: (i) to what 
extent the effects of the lattice persist as we go away from 
the half filling, and (ii) how the pairing crosses over to 
the pairing in the electron gas in the dilute regime. We 
shall show that (i) the pairing symmetries {(1^2 _y2^ dxy) 
reflecting lattice structure persist well away (n > 0.2) 



2 



from the half fihed band, and (ii) in the dilute (n < 0.2) 
regime the pairing symmetries are s and p in agreement 
with those for the electron gas. The key factors con- 
trolling the dominant symmetries are identified to be the 
structure (peak positions, intensities and widths) of the 
charge and spin susceptibilities (along with the size of 
the Fermi surface in the dilute regime). 

We first introduce an extended Hubbard model with 
long-range interactions (Fig. P), 

nn ^ 
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Here t(= 1) is the nearest-neighbor transfer, taken to 
be the unit of energy, and the off-site Coulomb repulsion 
extends up to the 12th neighbor for the square lat- 
tice as depicted in Fig^ where the magnitude of the nth- 
neighbor repulsion = {Vi = V) is taken to be in- 
versely proportional to the distance (e.g., V12 = V/V20). 
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FIG. 1: A square lattice with the long-range repulsion K 



We have previously extended^^ the FLEX formalism to 
the extended Hubbard model with the nearest-neighbor 
interaction, but it is difficult to extend this beyond the 
second neighbors. So we employ here a simplified FLEX, 
where full bubble diagrams and restricted ladder dia- 
grams (that include only the U term) are considered 
for the effective interaction in a self-consistent manner. 
So this is a kind of RPA with all the Green's functions 
dressed, which was adopted in Refs^ii^ for the extended 
Hubbard model with the nearest-neighbor interaction. 
The simplified FLEX still belongs to "conserving approx- 
imations" formulated by Baym and Kadanofl?^*^^. The 
FLEX neglects vertex corrections, but it has been known 
that the correction can lead to some problems in self- 
consistent calculations, especially in the dilute regime^ 
So here we opt for FLEX and assume that the vertex 
correction would not significantly affect the competition 
between different pairing symmetries (although Tc, not 
discussed here, may be affected). 

The spin (Xsp) and charge (Xch) susceptibilities are 

then Xsp{q) = i-u^{q) ' Xch{q) = i+[t/+2y(g)]x(g) ' 
where the irreducible susceptibility is given by x{q) = 
-{T/N)^^G{k ^ q)G{k). Here q = {q.Cn) with 
Cn = 2n7rT being the Matsubara frequencies for bosons 
and k = {k^uJn) with uJn = (2n — l)7rT for fermions. 



and V{q) = 2Vi(cosgaj + cosqy) + 4V2(cos cos g^y) + 
2F3(cos2ga^ + cos2qy) H . The self-energy is given as 

+ + 2UVik - k') + 2V^{k - /c')| Xch(fc - /c')| G{k').{2) 

The gap function and Tc are obtained with Eliashberg's 
equation, 

M{k) = E - k')G{k')G{-k')<i>{k'), (3) 

k' 

where A = 1 corresponds to T = Tc, and the pairing 
interactions rs(rt) for the singlet (triplet) channels are 

r.(g) = C/ + V{q) + ^f/'xsp - [if/' + 2UV{q) + 2^' (g)|xch (<?) ,(4) 
Vt{q) = V{q) - if/'xsp - [it/' + 2UV{q) + 2^' (g)|xch (?) • (5) 

Here we take A/" = 32 x 32 /c-point meshes and the Mat- 
subara frequencies uOn from — {2Nc — 1)ttT to (2Nc — 1)7tT 
with Nc = 16384, T = 0.01, U = 4.0. 

We first show the eigenvalue. A, of Eliashberg's equa- 
tion, along with the charge (Xch) and spin (Xsp) suscep- 
tibilities against the band filling n for V = 1.0 in Fig|21 
While A for T = 0.01 is still much smaller than unity, it is 
difficult to go down to lower temperatures, given a huge 
computational demand for the long-range model. So we 
assume that the dominant pairing symmetry is the one 
that has the largest A at T = 0.01 which is higher than 
Tc. This amounts to assuming that, while the position of 
boundary may shift as T ^ 0, the order in which dom- 
inant symmetry changes is not altered. We see that the 
dominant pairing is (1^2 _y2 d^y ^ p ^ s for the band 
filling is decreased from the half- filling. Figure |31 depicts 
the four gap functions in k space for typical values of n. 

To explore why the dominant gap symmetry changes 
in such a way, we have plotted the charge and spin sus- 
ceptibilities in k space in Fig. ^ Let us start with the 
region 0.8 > n > 0.3, where the symmetry changes from 
dx2_y2 to dxy as n is decreased. For n = 0.8 the peak 
positions in the spin susceptibility are situated around 
k = (tt, tt) as marked with an arrow in Fig. ^ which 
accounts for the dj;2_y2 pairing. For n = 0.5, the peaks 
shift to (0, ±7r), (±7r, 0), which should be why the dxy 
gap function, which changes sign across the pair hopping 
(0, ±7r), (±7r, 0), takes over. 

If we further decrease the band filling, we see from 
FigElthat Xch, which is originally much smaller than Xsp 
for n > 0.3, becomes comparable with Xsp- When the 
spin-fluctuation mediated interaction is dominant (as is 
the case with short-range repulsions), the singlet pairing 
interaction (eq^ has the Xsp-term whose coefficient is 
three times larger in magnitude than that for the triplet, 
F^ (eqEJ. So the singlet superconductivity is generally fa- 
vored when the interaction is short-ranged (i.e., the spin 
fluctuation is dominant). Conversely, if spin and charge 



3 



1 

0.6 
0.4 
0.2 



r=o.oi, u=4, v=\ 



P 



dxy 



triplet 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

n 



SDW 




0.1 0^ 0.3 0.4 0.5 0.6 0.7 O.S 0.9 

H 

FIG. 2: Top: Maximum eigenvalue, A, of Eliashberg's equa- 
tion in the triplet (solid line) and singlet (dotted) channels 
with the dominant symmetry indicated along with the bound- 
ary values of Vs on the top axis. The SDW phase is identi- 
fied from the divergence of the spin susceptibility. Bottom: 
Charge (solid line) and spin (dotted) susceptibilities against 
n. 



susceptibilities are comparable, the situation may be in- 
verted, which indeed occurs for the triplet p-wave around 
0.2 > n > 0.05. 
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FIG. 4: Spin susceptibility Xsp (top panels) and charge sus- 
ceptibility Xch (bottom) in k space for typical values of the 
band filling n. To make the peaks clearer, the color coding 
differs from frame to frame. 
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FIG. 3: Fermi surface (obtained as + ReS(fe) = /x; black 
line) and nodes in the dominant gap function (blue lines for 
singlet and red for triplet) for 0.04 < n < 0.8 with U = 4.0 
and V — 1.0. The arrow indicates the main scattering process 
mediated by spin fluctuations. 



FIG. 5: The maximum eigenvalue of Eliashberg's equation for 
triplet (solid line) and for singlet (dotted) against the range of 
the interaction (in units of the lattice constant) for n = 0.03 
and y = 1.0. 



In the most dilute regime (n < 0.05) an s-wave has 
the largest A. Since the ratio Xch/Xsp does not drasti- 
cally change in this region, the appearance of s requires 
an explanation. Here we can note another factor that 
affects the pairing. This is a very general observation 
that nodes in the gap function, even when they are nec- 
essary, act to lower Tc, since some of the pair scatterings 
around each node work against the paired state. In the 
present context, this occurs for the p wave for which the 
Fermi surface shrinks with decreasing n, so that the frac- 
tion of the phase space volume for the unfavorable region 
around the nodes increases. We identify this to be the 
reason why p gives way to s in the most dilute region. 
The competition between s and p when the range of in- 
teraction is varied is displayed in FigO where we see 
that s-wave tends to be favored as the interaction be- 
comes longer-ranged. We can understand this in terms 
of the pairing interaction, since the spin structure in this 
regime is peaked at = 0, which inhibits the singlet s 
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and assists the triplet and the peak rapidly decreases 
as the range is increased. 

Let us finally compare the present result with the sit- 
uation in the electron gas. Takada^ has obtained, with 
the Kukkonen-Overhauser method^^, a phase diagram as 
a function of the electronic-density parameter Ts for the 
electron gas, where p-wave (at = 3.3) and then 5- wave 
(at Ts = 8.6) appear as Ts is increased (i.e., as the gas be- 
comes more dilute). To quantify the connection between 
the electron gas and the present lattice model (Fig|21), we 
have to define for the latter. In the effective- mass sense 
we have = \/2/{a%kF), where a% = (47reof^^)/(m*e^) 
is the effective Bohr radius, and m* ~ h^/{2o?t) the 
effective mass with a being the lattice constant. Still, 
the definition can be a bit tricky, since the interaction is 
not truly 1/r. Here we eliminate e from the above and 
V = e^/(47reoa) to define = V/ {^/2akFt)^ where kp 
in the radius of the Fermi circle when the band is dilute, 
while for n close to half-filling we define /c^ as the largest 
radial distance from F as a measure of the Fermi surface. 
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FIG. 6: The critical value of Vs between s-wave and p-wave 
plotted against V/t for [/ = 4.0 at T = 0.01. 



We first examine how the boundary value of Vg between 
s- and p-waves changes with V in FigEl We can see that 
the boundary is almost flat, especially so for V/t < 1 
where the dispersion around small /ci? is close to those 
of the electron gas. So the result is consistent in that 
the phase is basically determined by Vg alone as in the 
electron gas. 

If we turn to the numerical values of at the bound- 
ary, we see that, although p- and 5- waves appear in the 
dilute concentration regime in both of the present and 
electron-gas^ models, the boundary values of Vg exhibit 
significant differences between the two. One factor for 
the discrepancy should be that the present model is a 
2D system, while Ref*^ considers 3D systems. In the 
plasmon-mechanism analysis of superconductivit}ii& the 
critical in 2D is seen to be reduced (about 1/3) from 
that for 3D. More precisely, the present model differs 
from the Coulomb gas even though the 12-th neighbor 
interaction is included, since a cut-off in the interaction 
should degrade Gauss's law (hence degrade the plasmon 
modes). Extending our model to 3D may cast a light on 
this. 
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